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ABSTRACT 

Context. We study the properties of the X-ray surface brightness profiles in a sample of galaxy clusters that were observed with 
Chandra and have emission detectable with a signal-to-noise ratio higher than 2 per radial bin at a radius beyond -R500 ~ 0.7 x -R200 ■ 
Aims. Our study aims to measure the slopes in both the X-ray surface brightness and gas density profiles in the outskirts of massive 
clusters. These constraints are compared with similar results obtained from observations and numerical simulations of the temperature 
and dark-matter density profiles with the intention of presenting a consistent picture of the outer regions of galaxy clusters. 
Methods. We extract the surface brightness profiles Sb(r) of 52 X-ray luminous galaxy clusters at z > 0.3 from X-ray exposures 
obtained with Chandra. These objects, which are of both high X-ray surface brightness and high redshift, allow us to use Chandra 
either in ACIS-I or even ACIS-S configuration to survey the cluster outskirts. We estimate R200 using both a /3— model that reproduces 
the surface brightness profiles and scaling relations from the literature. The two methods converge to comparable values. We determine 
the radius, Rs2n, at which the signal-to-noise ratio is higher than 2, and select the objects in the sample that satisfy the criterion 
Rs2N / R200 > 0.7. For the eleven selected objects, we model by a power-law function the behaviour of Sb{r) to estimate the slope 
at several characteristic radii expressed as a fraction of -R200 • 

Results. We measure a consistent steepening of the Sb{r) profile moving outward from 0.4_R200, where an average slope of —3.6 
(cr — 0.8) is estimated. At -R200, we evaluate a slope of —4.3 (a — 0.9) that implies a slope in the gas density profile of ~ —2.6 and 
a predicted mean value of the surface brightness in the 0.5 — 2 keV band of 2 x 10~^^ erg s~^ cm~^ deg~^. 

Conclusions. Combined with estimates of the outer slope of the gas temperature profile and expectations about the dark matter 
distribution, these measurements lie well within the physically allowed regions, allowing us to describe properly how X-ray luminous 
clusters behave out to the virial radius. 

Key words, galaxies: cluster: general - intergalactic medium - X-ray: galaxies - cosmology: observations - dark matter. 



1 . Introduction 

Galaxy clusters form by the hierarchical accretion of cosmic 
matter The end products of this process are virialized struc- 
tures that, in the X-ray band, exhibit similar radial profiles of 
surface brightness (e.g. Vikhlinin et al. 1999, Neumann 2005), 
plasma temperature (e.g. Allen et al. 2001, Vikhlinin et al. 2005) 
and dark matter distribution (e.g. Pointecouteau, Arnaud & Pratt 
2005). These measurements have been improved due to the arc- 
sec resolution and large collecting area of the X-ray satellites, 
such as Chandra and XMM-Newton, but still remain difficult be- 
cause of the high signal-to-noise ratio required. On the other 
hand, the X-ray surface brightness is a far easier quantity to 
observe and define, which is rich in physical information be- 
ing proportional to the emission measure of the emitting source. 
Recent work focused on a few local bright objects with ROSAT 
PSPC observations, which have low instrumental background 
and large field of view, to recover and characterise the X-ray 
surface brightness profile over a significant fraction of the virial 
radius (Vikhlinin et al. 1999, Neumann 2005). 

In this work, we study the X-ray emission of a sample of 
galaxy clusters in the redshift range 0.3 — 1.3 observed with the 
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arcsec resolution of Chandra (see Balestra et al. 2007), with the 
purpose of mapping, even at high redshift, a region out to i?2oo 
and place reasonable constraints on the slope of the gas density 
in the outskirts of hot (kT > 3 keV) galaxy clusters. 

We assume a Hubble constant of 70 hro km s^^ Mpc^^ in 
a flat universe with fl^ equal to 0.3. All quoted errors are pre- 
sented as the 1 sigma confidence level, unless otherwise stated. 

2. Preparation of thie dataset 

To study the outskirts of X-ray emitting galaxy clusters, possi- 
bly out to the virial region, we require bright objects located at 
an appropriate redshift to allow the Chandra field-of-view to en- 
compass the interesting area. As shown in Fig. [T] (panel on the 
left) for a typical massive object, we should expect to resolve 
the region within an overdensity of 200, estimated with respect 
to the critical density at the cluster redshift, using instruments 
with a field-of-view of radius larger than 8 arcmin above red- 
shift 0.2. In the same figure, we plot the expected dependence 
of the radius (as a function of i?2oo) on both the cluster over- 
density (central panel) and the gas temperature and redshift for 
a given surface brightness (panel on the right). These plots in- 
dicate that any overdensity is mapped within a fixed fraction of 
R2Q0 regardless of the cluster temperature and/or mass concen- 
tration (e.g. Ac = 500 is reached at about 0.65i?2oo for any T 
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Fig. 1. (Left) Expected i?200 for an object with 4 (dashed line) and 8 keV (solid line) and radius corresponding to 4 and 8 (thickest 
line) arcmin for the assumed ACDM cosmology as function of the observed redshift. Two estimates, almost indistinguishable, of 
i?200 provided, one obtained with a /3— model with core radius of 0.1 Mpc and (3 — 0.6 (thickest lines), the other from the 
best-fit relation in Arnaud et al. (2005): ft.^i?200 = 1714(±30)(r/5keV)°-5°=^" °5. (Center) Fraction of i?2oo mapped within a 
given overdensity Ac at redshift for a cluster with a gas temperature of 4 (diamonds and dashed line) and 8 keV (solid line). The 
region at p/pc = 500 is enclosed within a sphere with radius « 0.65 x i?2oo- The dotted line indicates the mapped fraction of i?200 
assuming a Navarro-Frenk- White profile (e.g. Navarro et al. 2004) with concentration of 6, typical for the massive objects under 
consideration. (Right) Regions mapped (as fraction of i?20o) for fixed surface brightness at different redshifts and gas temperatures, 
Sb ^ T°-^(l + z)^*. As reference, a value of Sb = 1 at R200, z — 0, kT — 8 keV and a /?— model with core radius of 0.1 Mpc and 
/? = 0.6 are adopted. 

Table 1. Sample of galaxy clusters, which satisfies the criterion Rs2N / R20Q > 0.7. The columns show: the name of the cluster, 
the observation IDs, in which ACIS configuration has been observed, the exposure time, the adopted X-ray center, the Galactic 
absorption from Dickey & Lockman (1990) in correspondence of the X-ray center, the redshift, the gas temperature, the mean local 
background, B, in counts observed in the 0.5 — 5 keV band, and the average distance Db of the regions used to estimate the local 
background from the X-ray center Note: MACSJ0744.9H-3927 is the merger of three exposures with IDs 3197, 3585, and 6111. 



cluster 


Obs ID 


ACIS 


^cxp 


RA 


Dec 




z 


kT 


B(xl0"^) 


Db 








ks 


h m s 


degree 


lO^" cm-3 




keV 


cts/s/arcmin^ 


R200 


MS0015.9-F1609 


520 


I 


67.4 


18 33.8 


+ 16 26 12 


4.1 


0.541 


9.59 ± 0.43 


1.61 ±0.03 


1.35 


MACSJ2228.5+2036 


3285 


I 


19.9 


22 28 33.9 


+20 37 15 


4.6 


0.412 


8.25 ± 0.59 


2.00 ±0.07 


0.92 


MACS0744.9+3927 


Merged 


I 


89.0 


7 44 52.8 


+39 27 26 


5.7 


0.686 


9.58 ±0.67 


1.67 ±0.04 


1.64 


MACSJ0417.5-1154 


3270 


I 


11.8 


4 17 34.6 


-11 54 33 


3.9 


0.440 


10.84 ± 0.98 


2.33 ±0.10 


0.84 


RXJ1701.3+6414 


547 


I 


49.5 


17 01 23.8 


+64 14 1 1 


2.6 


0.453 


4.36 ± 0.28 


2.25 ±0.05 


1.18 


MACSJ1720.2+3536 


3280 


I 


20.8 


17 20 16.8 


+35 36 26 


3.4 


0.391 


6.46 ± 0.33 


1.87 ±0.09 


0.91 


RXJ1416.4+4446 


541 


I 


28.5 


14 16 27.8 


+44 46 45 


1.2 


0.400 


3.43 ± 0.20 


2.22 ±0.06 


1.38 


MACSJ2129.4-0741 


3199 


I 


17.6 


21 29 26.6 


-7 41 29 


4.8 


0.570 


9.17 ±0.90 


1.88 ±0.09 


1.12 


MACSJ1621. 3+3810 


3254 


I 


9.8 


16 21 24.9 


+38 10 08 


1.1 


0.465 


6.62 ± 0.74 


1.58 ±0.10 


1.21 


MS045 1.6-0305 


902 


S 


42.6 


4 54 11.3 


-3 00 56 


5.0 


0.540 


9.11 ± 0.45 


3.98 ± 0.07 


0.76 


MACSJ1206.2-0847 


3277 


I 


234 


12 06 12.1 


-8 48 02 


3.7 


0.440 


11.98 ±0.85 


2.07 ±0.06 


0.78 



of interest, and with variations of 2 per cent with a concentration 
parameter in the range 4 — 8 for a given Navarro-Frenk- White 
mass profile; e.g. Navarro et al. 2004), and the surface brightness 
that a galaxy cluster is expected to emit for a given T and z with 
respect to the value measured at z = in an object with T = 5 
keV (e.g. the same Sh value is expected in a 10 keV system at 
2 ~ 0.3). 

In the present work, we consider a sample of hot (Tgas > 3 
keV), high-redshift (0.3 < z < 1.3) galaxy clusters described 
in Balestra et al. (2007). We refer to that work for details of 
the data reduction. We recall that the spectral analysis was per- 
formed by extracting the spectrum from a circular region of ra- 
dius ^ 0.15 — 0.3i?2oo defined in order to maximize the signal- 



to-noise ratio in each cluster These regions contain the core 
emission, which has not been removed due to the low count 
statistics for the high redshift systems under consideration. On 
the other hand, results from observations and numerical simula- 
tions (e.g. Santos et al. 2008, Ettori & Brighenti 2008) indicate 
a lower incidence of cooling cores in clusters at higher redshift, 
suggesting that our overall estimates of Tgas should not be signif- 
icantly biased towards low values. Moreover, the measurements 
of the gas temperature were used in our analysis only to infer 
a physical radius, which depends on the square root of Tgas- 
Therefore, any error affecting our temperature estimates prop- 
agates only half of its value as a relative error on the physical 
radius of interest. 
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We prepare the exposure-corrected images in the energy 
band 0.5 — 5 keV as described in Tozzi et al. (2003) and Ettori et 
al. (2004). The exposure maps for each cluster were computed 
by combining different instrumental maps at given energy, using 
weights defined from a thermal spectrum with the best-fit param- 
eters obtained from our spectral analysis. The variations in the 
exposure map were not expected to be significant in the adopted 
energy bandQ. We masked any detectable point sources with cir- 
cular regions of radii large enough to include all the emission 
and, in any case, larger than 2 arcsec. We preferred to remove 
the point sources manually to avoid any possible confusion with 
detections by automatic algorithms (i.e. wavedetect) in proxim- 
ity or within the extended emission of the clusters. 

The center of the cluster was defined to correspond to the 
centroid of the raw image estimated in a box of width ^10 — 20 
arcsec around the maximum of the image itself that had been 
smoothed with a moving average in a square box of width equal 
to 3 — 5 arcsec. 

The surface brightness profiles were extracted by requiring 
a fixed number of a minimum 50 counts per bin in the inner 
10 radial bins and then collecting the counts within annulii for 
which the outer radius was increased at each step by a factor of 
1.2. 

A local background, B, was defined for each exposure by 
considering a region far from the X-ray center that covered a 
significant portion of the exposed CCD with negligible cluster 
emission. The estimated average value of B, its relative error, 
and the mean distance from the X-ray center of the region anal- 
ysed to determine the background are quoted in Table [1] 

We define the "signal-to-noise" ratio, S2N, to be the ra- 
tio of the observed surface brightness value in each radial bin, 
Sb{r), after subtraction of the estimated background, B, to 
the Poissonian error in the evaluated surface brightness, eb(r), 
summed in quadrature with the error in the background, e^: 
S2N{r) = [St{r) - B] / y/eb{r)^ + ej^. The outer radius at 
which the signal-to-noise ratio remained above 2 was defined 
to be the limit of the extension of the detectable X-ray emission. 

We note that the local background is about a factor of 2 
higher in the ACIS-S field than in AClS-1. We observe no signif- 
icant correlation of B with the Galactic absorption (< 6 x 10^° 
cm~^ in the selected objects; a random deviation from the null 
value of no-correlation is expected with probability of 0.19), 
indicating that we were dominated by cosmic and instrumen- 
tal background. A negligible (with probability of 0.69) correla- 
tion is also noticed between B and the mean distance Db (see 
Table [TJ of the region from which the local background is ex- 
tracted. A slightly positive trend (significance of 1.6 x 10~^) 
was instead present between the estimated uncertainty es and 
the exposure time of the observation, as expected. 

To scale the radial quantity with respect to a physical radius, 
we estimated i?2oo by using the single temperature measure- 
ment (see Table (TJ and the best-fit description of Sb{r) with a 
/3— model (the best-fit parameters are shown in Table|2] obtained 
after fitting the background-subtracted profile over all positive 
values above 40 kpc) following the equation in Ettori (2000): 

' See, for example. Fig. 2a in 

|http://cxc.harvard.edu/ciao3.4/download/doc/expmap-intro.ps| 



where = [^^(1 + z)^ + 1 - f^J " ■^ = 200, a = 

[1.5/3(7 ^ 1) + 1] and the polytropic index 7 is fixed to be 
equal to 1 . 

For comparison, we also estimated i?200 by adopting 
the scaling relation with Tgas obtained for massive clus- 
ters in Arnaud et al. (2005; si milar results were pre- 
sented in VikhUnin et al. 2005): Ihoo = 1714(±30) x 
(rgas/5keV)°-5°±'^ 05/i7oii;-i kpc. The value of was es- 
timated and quoted in Table |2] No appreciable differences were 
present between the different estimates of i?200i such that the 
mean and standard deviation of the ratio of the two values were 
1.02 and 0.06, respectively. 

3. The surface brightness at r > R500 

We investigated the X-ray surface-brightness profiles of massive 
clusters at r > i?5oo ~ 0.7i?2oo (see Fig. [TJ, selecting the 11 
objects with Rs2N / ^200 > 0.7. The properties of these clusters 
are shown in Tables [T] and |2l Examples of the analysed dataset 
are shown in Fig.|2] 

We performed a linear least-squares fit between the logarith- 
mic values of the radial bins and the background-subtracted X- 
ray surface brightness, taking into account the data errors in both 
coordinates with the routine FITEXY (Press et al. 1992; as im- 
plemented in IDL). 

The distribution of the error-weighted mean slopes above 
a fixed fraction of i?2oo and within Rs2N is shown in Fig. |3] 
and quoted in Table[3] Overall, the error-weighted mean slope is 
—2.91 (with a standard deviation in the distribution of 0.46) at 
r > 0.2i?2oo and -3.59(0.75) at r > 0.4i?2oo- For the only 3 
objects for which a fit between 0.5i?2oo and Rs2N was possible, 
we measured a further steepening of the profiles, with a mean 
slope of —4.43 and a standard deviation of 0.83. 

We also fitted linearly the derivative (from numerical dif- 
ferentiation using three-point Lagrangian interpolation as de- 
scribed in Hildebrand 1987 and implemented in the IDL func- 
tion DERIV) of the logarithm Sb{r) over the radial range 
0.1i?2oo ~ Rs2N, excluding in this way the influence of the 
core emission (see values sq and si in Table |3]l. The average 
(and standard deviation a) values of the extrapolated slopes are 
then -3.15(0.46), -3.86(0.70), and -4.31(0.87) at 0.4i?2oo, 
0.7i?2oo and i?200, respectively. 

These values are comparable to these obtained in previous 
analyses. Vikhlinin et al. (1999) found that a /3— model with f3 = 
0.65 — 0.85 described the surface brightness profiles in the range 
0.3 — I-R180 of 39 massive local galaxy clusters observed with 
ROSAT PSPC. These values translate in a range of the estimate 
of the power-law slope(-2.9/-4.1,from56(r) ~ ,-2(0.5-3/3) ^ 
j,i-6/3-) jjj^j. includes our estimates. 

Neumann (2005) found that the stacked profiles of a few 
massive nearby systems located in regions of low (< 6 x 10^" 
cm^'^) Galactic absorption observed by ROSAT PSPC still pro- 
vided values of (3 around 0.8 at R20Q, with a power-law slope 
that increased from —3, when the fit was done over the radial 
range [0.1, l]i?2oo, to -5.7^^;^ over [0.7, 1.2]i?2oo. 

These observational results were supported by hydrodynami- 
cal simulations of X-ray emitting galaxy clusters performed with 
the Tree-nSPH code GADGET-2 from Roncarelli et al. (2006). In 
the most massive systems, they measured a steepening of 5f,(r), 
independently from the physics adopted in treating the baryonic 
mass, with a slope of —4, —4.5, and —5.2 when estimated for the 
radial range 0.3 - 1.2i?200, 0.7 - 1.2i?2oo, and 1.2 - 2.7i?2oo, 
respectively. In particular, we note the good agreement between 
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Fig. 2. In the first column, we plot the surface brightness profiles of the objects under examination with the fitted background 
(horizontal dotted line) and the radius i?2oo (vertical dashed line). The panels in the second column show the signal-to-noise 
profiles, evaluated as 5*2 = (5*6 — B)/e, where the error e is the sum in quadrature of the Poissonian error in the radial counts 
and the uncertainties in the fitted background, B. They are plotted as a function of r/i?2oo- In the third column, we plot the best-fit 
values of the slope of the surface brightness profiles as a function of r/i?2oo- These values are estimated over 6 radial bins. The thick 
horizontal solid line indicates the slope evaluated between 0.4 x i?200 and Rs2N with a minimum of 3 radial bins. The dashed line 
indicates the best-fit of dln{Sb) / dln{r / R200) with the functional form sq + si ln(r/i?20o) over the radial range 0.1 x R200 — Rs2N , 
with the best-fit parameters quoted in Table [?] 



the slope of the simulated surface brightness profile of the rep- 
resentative massive cluster in the radial bin 0.7 — 1.2_R200 (see 
values of 6^ in Table 4 of Roncarelli et al. (2006), ranging be- 
tween —4.29 and —4.54) and our mean extrapolated value at 
i?200 of -4.31. 



In Table [3] we quote the estimated surface brightness value 
at i?200 in CGS units. They were obtained by extrapolating the 
measured St, at ^ 0.7i?2oo out to i?200 by using the best- fit val- 
ues quoted in the same Table. The conversion from the count 
rate to the flux was completed by adopting an absorbed ther- 
mal component at the cluster redshift with an assumed metal- 
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Table 2. Best-fit results of the X-ray surface brightness profiles modelled with a j3- 
in Arnaud et al. (2005); i?2oo from Eq.lT] Rs2N and are quoted in unit of -R2oo- 



model. i?200 is estimated from the scaUng laws 



cluster 



-R200 
kpc 



R200 
kpc 



RS2N 
R2OO 



r-c 
R200 



MS0015.9+1609 


1770 


1948 


± 


43 





871 





137 


± 





005 





742 


± 





oil 


MACSJ2228.5+2036 


1771 


1743 


± 


62 





837 





093 


± 





007 





587 


± 





013 


MACS0744.9+3927 


1622 


1566 


± 


56 





832 





045 


± 





003 





561 


± 





005 


MACS J0417.5- 1154 


1998 


1936 


± 


91 





783 





068 


± 





006 





567 


± 





010 


RXJ1701.3+6414 


1258 


1204 


± 


44 





731 





087 


± 





013 





555 


± 





020 


MACSJ1720.2+3536 


1587 


1617 


± 


42 





726 





056 


± 





004 





626 


± 





012 


RXJ1416.4+4446 


1150 


1133 
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38 





722 





080 


± 





010 





587 


± 





020 


MACSJ2129.4-0741 


1701 


1684 
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83 





712 





069 


± 





006 





592 
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014 
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1538 


1516 


± 


88 





711 





048 


± 





007 





585 


± 





019 


MS045 1.6-0305 


1727 


1995 


± 


52 





705 





118 


± 





005 





814 


± 





015 
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Fig. 3. Error-weighted mean and standard deviation of the distri- 
bution of the values of the slope of the surface brightness in the 
radial range r — Rs2N plotted as a function of the inner radius r. 
The number of objects with at least 3 points considered for this 
estimate is indicated. 



licity and temperature equal to (i) the ones measured and (ii) 
one third of the measured values. We predict an average sur- 
face brightness of about 2 x 10"^^ erg s~^ cm~^ deg~^ in the 
0.5 — 2 keV band, corresponding to about half of the value asso- 
ciated with the total cosmic background observed (e.g. Hickox 
& Markevitch 2006) and consistent with the estimates measured 
in the hydrodynamical simulations of massive galaxy clusters 
(Roncarelli et al. 2006). Converting this surface brightness to 
counts in the 0.5 — 5 keV band, we predict that at i?200 the X- 
ray emission from a cluster is responsible of 5 per cent of the 
total counts observable with Chandra. 



4. Physical properties of the outer cluster 
regions 

We combine the measured slope of the surface brightness in the 
cluster outskirts with the predicted slopes of the gas temperature 
and dark matter profiles. We have demonstrated that the slope 
depends on the radius at which it is measured. Thus, we focus 
hereafter on the properties at R200, by representing with power- 
laws the radial dependence of the gas density, Ugas, temperature, 
Tgas, and dark matter, pdm, profiles: 



PUM OC r" 



(2) 



The surface brightness, St, is the integral along the line of 
sight of the square of the gas density and can be written, assum- 
ing spherical geometry, as 



.2a +1 +0.26 



(3) 



where we consider the weak dependence (T oc r^-^^) on the gas 
temperature of the cooling function, A(r), due to the limited 
energy band AE = [0.5, 5] keV in which the images were con- 
sidered and the gas temperatures under consideration (e.g. Ettori 
2000). 

Hydrostatic polytropic gas can be described by the relation 
Tgas oc nJ^J-, with a polytropic index 7 that relates the ratio of 
specific heats in a adiabatic gas and equals 5/3 for a monoatomic 
ideal gas. If the intracluster gas is not entirely adiabatic, the poly- 
tropic index can be used as a fitting parameter to describe the 
relative behaviour of the gas temperature and density profile. In 
particular, if the thermal conduction is an efficient process, the 
gas tends to be isothermal and 7 = 1, whereas the entropy per 
atom will be constant if the thermal conduction is slow and the 
gas is well mixed. Moreover, when 7 > 5/3, the gas becomes 
convectively unstable and mixing occurs within several sound 
crossing times, which are in general quite short compared to the 
overall age of the cluster Therefore, this "polytropic condition" 
requires that 1 < 7 < 5/3 and translates into the request that 



b 2 

0< - < 
a 3 



(4) 
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Table 3. Best-fit values of the slope of the surface brightness profile in the radial range [0.2 x -R200, Rs2n] and [0.4 x R2Q0, Rs2n]- 
The best-fit results of the first derivative of Sb{r) over the entire radial range are indicated as d\n{Sb) / d\n{r / R200) = sq + 
si ln(r/i?2oo) and shown in Fig. |2l The predicted slopes from the best-fit results at 0.4, 0.7, and 1 x i?200 are quoted. In the last 
two columns, we present the expected surface brightness in the 0.5 — 2 keV band and the fraction /c2oo of the clusters counts with 
respect to the total count rate measured in the 0.5 — 5 keV at i?2oo- 



cluster 


slope 

0.2_R200 — RS2N 0.4_R200 — -RS2JV 


(so, Sl) 


0.4i?2oo 


0.7i?2oo 


-^200 


*5'b(-^20o) 

10~^^ erg/s/cm^/deg^ 


/C200 


MS0015.9-F1609 


-3.30 ±0.28 


-3.89 ±0.96 


(—4.41^ -1.33) 


-3.19 


-3.93 


-4.41 


(1.7, 1.7) 


0.051 


MACSJ2228.5-^2036 


-2.74 ±0.29 


-4.26 ± 1.23 


(-4.92^-1.75) 


-3.32 


-4.30 


-4.92 


(2.6^2.7) 


0.057 


MACS0744.9+3927 


-2.73 ±0.31 


-3.23 ±0.89 


(-3.73, -0.88) 


-2.92 


-3.41 


-3.73 


(2.6,2.7) 


0.074 


MACSJ0417.5-1154 


-3.38 ±0.42 


-4.53 ± 1.27 


(-5.60,-2.02) 


-3.75 


-4.88 


-5.60 


(1.6, 1.6) 


0.034 


RXJ1701.3+6414 


-2.39 ±0.27 


-2.59 ± 1.17 


(-4.24,-1.60) 


-2.78 


-3.67 


-4.24 


(2.4,2.6) 


0.051 


MACSJ1720.2+3536 


-3.19 ±0.36 


-2.40 ± 1.32 


(-3.57,-0.77) 


-2.87 


-3.30 


-3.57 


(2.5,2.6) 


0.059 


RXJ1416.4+4446 


-2.74 ±0.37 


-3.63 ±1.72 


(-3.75, -0.94) 


-2.89 


-3.41 


-3.75 


(1.8,1.9) 


0.042 


MACSJ2129.4-0741 


-2.68 ±0.33 


-4.16 ± 1.74 


(-3.92,-1.13) 


-2.88 


-3.51 


-3.92 


(2.4,2.5) 


0.060 


MACSJ1621. 3+3810 


-2.64 ±0.37 


-3.53 ± 1.76 


(-2.99, -0.45) 


-2.58 


-2.83 


-2.99 


(3.0,3.2) 


0.087 


MS045 1.6-0305 


-3.94 ±0.45 


-4.77 ± 2.04 


(-5.85,-1.91) 


-4.11 


-5.17 


-5.85 


(0.3,0.3) 


0.006 


MACSJ1206.2-0847 


-3.24 ±0.37 


-3.67 ± 1.64 


(-4.43,-1.16) 


-3.37 


-4.02 


-4.43 


(1.4,1.5) 


0.032 


mean 


-2.91 


-3.59 


(-4.31,-1.27) 


-3.15 


-3.86 


-4.31 


(2.0,2.1) 


0.050 


a 


0.46 


0.75 


(0.87,0.50) 


0.46 


0.70 


0.87 


(0.8,0.8) 


0.022 



Further constraints are provided by the integrated masses. 
The gas mass estimate. A/gas ~ / UgasdV oc r^^^, requires 
that 

a ± 3 > 0. (5) 

The total gravitating mass can be expressed both as an integral 
of the dark matter profile over the cluster volume 



A4 



= / PDudV oc 



„c+3 



(6) 



and with the assumption of hydrostatic equilibrium between dark 
matter and spherically symmetric intracluster gas 



hyd 



-rT 



din 'I 



^gas 



(Jin?' 



d In Tga, 
dlnr 



oc — (a±6) r 



b+l 



(7) 



The request for a gravitational mass increasing radially implies 
that c + 3 > and a + b < 0. Moreover, by equating the two 
definitions of the gravitational mass, we find that the exponents 
of the radial dependence must to satisfy the condition c + 3 = 
6+1, which directly relates the slope of the dark matter profile 
to that of the gas temperature profile. 

Overall, these conditions can be summarized by the follow- 
ing inequalities among the indices of the assumed power-laws: 



(8) 



We show in Fig.|4]the allowed regions in the a — b — c plane. 

It is worth noting that the permitted values shown in these 
plots also satisfy the requests that the sound-crossing time, 

tsc oc TgJs^^r « ri~^/2, the cooHng time, tcooi oc TglsU'^ « 

and the equipartition time for Coulomb collisions, teq oc 

Tglsn~^g w increase with radius in the cluster out- 

skirts. 




Some observational constraints are available and can be lo- 
cated in these graphs: (1) the mean value of the slope in the sur- 
face brightness at i?2oo is —4.31 (a = 0.87), equal to 2a ± 1 + 
0.2b as inferred by Eq. [3] and implies that a — —2.62 ± 0.43; 
(2) the dark matter profile behaves as 7--2.5±o.2 discussed and 
shown in Fig. 3 of Navarro et al. (2004); (3) a polytropic index 
7 « 1.24 is measured in the outskirts of temperature profiles 
(e.g. Markevitch et al. 1998, De Grandi & Molendi 2002 out to 
about 0.6i?i8o ~ 0.63i?2oo; a value of 1.42 ± 0.03 is estimated 
in one of the first spectral determination of the temperature pro- 
file at i?2oo, obtained for the cluster PKS0745-191 with Suzaku 
in George et al. 2008). These further constraints are also shown 
in Fig. i 



The expected slope of the dark matter profile, combined with 
the observed slope in surface brightness, places constraints on 
the predicted behaviour of the temperature profile (Eq.|8]l and its 
polytropic index, 7 = 1 + fe/a. For the values quoted above, we 
expect that b « -0.5 ± 0.2 and 7 « 1.19 ± 0.08. These values 
agree well with the constraints obtained from observed data-sets 
(e.g. Markevitch et al. 1998, VikhUnin et al. 2005, Leccardi & 
Molendi 2008, Reiprich et al. 2008, George et al. 2008), hydro- 
dynamical simulated objects (e.g. Loken et al. 2002, Roncarelli 
et al. 2006) and analytical model of the ICM (e.g. Ostriker, Bode 
& Babul 2005). 



Finally, the entropy profiles, K{r) — rgas(?')/"-gas, both 
from smooth and hierarchical accretion models (e.g. Tozzi & 
Norman 2001, Voit 2005) are predicted to increase with radius 
as K oc r^-^. This shape of the entropy profile outside the core, 
but well within the virial radius, is also observed in high qual- 
ity cluster exposures with the XMM-Newton satellite (Pratt & 
Arnaud 2003). By using our power-law expressions, we find 
that K{r) oc 7-''-2a/3^ \N'\\h a predicted value of the slope of 
1.25 ± 0.35 and 1.12 for the simulated dark matter distribution 
and the observed polytropic index, respectively. 
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1.0 1.5 2.0 2.5 3.0 

-a 



Fig. 4. Constraints on the outer slopes of the gas density (index 
a), temperature (index b) and dark matter (index c — b ~ 2) pro- 
files. The shaded area shows the physically allowed regions for 
the index values as discussed in the text. The hatched regions 
indicate the la constraints from Navarro et al. (2004; horizon- 
tal regions) and from the present work including the dependence 
upon the cooling function (vertical area). The horizontal green 
solid line shows the fiducial slope of the temperature profile from 
Vikhlinin et al. (2005). The vertical solid (and dashed) lines indi- 
cate the limits on the slope of the gas density profile from hydro- 
dynamical simulations (tighter constraints in orange; Roncarelli 
et al. 2006) and from the present observational analysis (in red). 
The diagonal dashed line is the locus of the predicted behaviour 
of the entropy profile, b — 2a/ 3 = 1.1. The dot-dot-dot-dashed 
line satisfies the condition that the gas fraction is constant in the 
cluster outskirts. 



5. Summary and Conclusions 

We have selected 1 1 massive galaxy clusters in the redshift range 
0.3 < z < 0.7 observed with Chandra, which is convenient to 
survey the X-ray emission out to a significant fraction of the 
virial radius and, in all the cases, to radii beyond i?5oo ~ 0.65 x 

-^200- 

By fitting a single power-law function to the X-ray surface 
brightness profile 5f,(r) in different radial ranges, we have de- 
tected a consistent steepening in Sb (r) moving outward, with a 
mean slope of —2.9 ±0.5 when the radial range 0.2i?2oo ^ Rs2N 
is considered, and —3.6 (r.m.s. 0.8) when the lower limit of the 
radial range is 0.4i?2oo- The mean slope estimated from the lin- 
ear fit to the derivative of Sb{r) is —3.2, —3.9, and —4.3 at 
0.4i?2oo, 0.7i?2oo, and i?200i respectively. These values, cor- 
rected by the weak dependence of the cooling function on gas 
temperature in the considered energy band [0.5, 5] keV, imply 
a slope in the gas density profile of -2.04 ± 0.23 (/3 « 0.68), 
-2.39 ± 0.35 (13 « 0.80) and -2.62 ± 0.43 (/? « 0.87) at 
0.4i?2oo, 0.7i?2oo and i?200, respectively. 

Combined with previous estimates of the outer slope of the 
gas temperature profile (e.g. Vikhlinin et al. 2005, from an ob- 



servational point of view; Roncarelli et al. 2006, for hydrody- 
namically simulated clusters) and the expectations for the dark 
matter distribution (e.g. Navarro et al. 2004), and modelling of 
the gas density, temperature and dark matter profiles with single 
power laws (rtgas oc r^jTgas oc r'',pDM « r"^), we have de- 
fined a physically-allowed region in the a — b — c plane, with a 
predicted behaviour at i?200 described by power-law indices of 
a « -2.6, b « -0.5, and c « -2.5. 
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